
Current Weather and Winds Aloft are shown alongside a pilot-selected route. Wind velocity at the 
pilot-selected altitude is depicted graphically with black arrows. The current weather is shown using 
symbolic and textual representations. 


human-centered methodology oriented 
towards providing the weather informa- 
tion (1) that the pilot needs and/or 
wants, (2) at the appropriate time, and 
(3) in the appropriate format. 

AWE can be characterized as a con- 
text-aware, domain-and-task knowledge- 
able, personalized, adaptive assistant. 


AWE automatically monitors weather re- 
ports for the pilot’s flight route and 
warns the pilot of any weather condi- 
tions outside the limits of acceptable 
weather conditions that the pilot has 
specified in advance. AWE provides tex- 
tual and/or graphical representations 
of important weather elements overlaid 


on a navigation map (see figure). The 
representations depict current and fore- 
cast conditions in an easy-to-interpret 
manner and are geographically posi- 
tioned next to each applicable airport 
to enable the pilot to visualize condi- 
tions along the route. In addition to au- 
tomatic warnings, the system enables 
the pilot to verbally request (via the 
speech-based user interface) weather 
and airport information. 

AWE is context-aware in the follow- 
ing sense: From the location of the air- 
craft (as determined by a Global Posi- 
tioning System receiver) and the route 
as specified by the pilot, AWE deter- 
mines the phase of flight. In determin- 
ing the timing of warnings and the 
manner in which warnings are issued, 
AWE takes account of the phase of 
flight, the pilot’s definition of accept- 
able weather conditions, and the pilot’s 
preferences for automatic notification. 
By noting the pilot’s verbal requests for 
information during the various phases 
of flight, the system learns to provide 
the information, without explicit re- 
quests, at the corresponding times on 
subsequent flights under similar condi- 
tions. 

This work was done by Lilly Spirkovska of 
Ames Research Center and Suresh K. 
Lodha of the University of California. Fur- 
ther information is contained in a TSP (see 
page 1). 

Inquiries concerning rights for the com- 
mercial use of this invention should be ad- 
dressed to the Patent Counsel, Ames Re- 
search Center, (650) 604-5104. Refer to 
ARC-14970-1. 


t? Fast Lossless Compression of Multispectral-Image Data 

A low-complexity adaptive-filtering algorithm is used. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


An algorithm that effects fast lossless 
compression of multispectral-image 
data is based on low-complexity, proven 
adaptive-filtering algorithms. This algo- 
rithm is intended for use in compress- 
ing multispectral-image data aboard 
spacecraft for transmission to Earth sta- 
tions. Variants of this algorithm could 
be useful for lossless compression of 
three-dimensional medical imagery 
and, perhaps, for compressing image 
data in general. 

The main adaptive-filtering algorithm 
on which the present algorithm is based 
is the sign algorithm (also known as the 


sign-error algorithm and as the binary 
reinforcement algorithm). The sign al- 
gorithm is related to the least-mean- 
square (LMS) algorithm. Both algo- 
rithms are briefly described in the 
following two paragraphs. 

Consider a sequence of image data 
(or any other data) that one seeks to 
compress. The sequence is specified in 
terms of a sequentially increasing index 
( k ) and the value (d k ) of the Ath sample. 
An estimated value of the Mi sample, d h 
is calculated by the equation 
d k = w*u*, 

where w* is a filter-weight vector at 


index k and u* is an input vector that 
can be defined in any of a number of 
different ways, depending on the spe- 
cific application. Once the estimate d k 
has been calculated, the error between 
the estimate and the exact value is calcu- 
lated as 

ek = d k - d k . 

When the LMS algorithm or sign algo- 
rithm is used as part of a predictive com- 
pression scheme, the sequence of e h val- 
ues is encoded in the compressed 
bitstream. 

The error value is also used to update 
the filter weights in either of two ways, de- 
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pending on which algorithm is in use. In 
the LMS algorithm, the update equation is 
w* + i = w* - pu k e k . 

In the sign algorithm, the update equation is 
w* + i = w*- /iujsgn(ej). 

In both update equations, R is a pos- 
itive, scalar step-size parameter that 
controls the trade-off between conver- 
gence speed and average steady-state 
error. A smaller value of R results in 
better steady-state performance but 
slower convergence. In some variants 
of these algorithms, the value of R 
changes over time. 

In the present algorithm, the index k 
is taken as an abstract representation of 
three indices ( x,y,z ) that are the coordi- 
nates of the sample in the multispectral 
dataset. Specifically, x and y are the spa- 
tial coordinates and z denotes the spec- 
tral band. The signal level (equivalently, 
the sample value) for that location is 
represented by d k = s(x,y,z). 

For purposes of compression, an image 
represented by a stream of data to be com- 
pressed is partitioned spatially into conve- 


niently sized, fixed regions. The data are 
compressed in the order in which they are 
received, maintaining separate statistics for 
each band and switching among the bands 
as necessary. The data from each region are 
compressed independently of those from 
other regions. Performing independent 
compression calculations for each region 
limits the adverse effect of loss of data. 

The input vector u,, chosen for this al- 
gorithm contains values from a six-sam- 
ple prediction neighborhood of a sam- 
ple of interest: three values from 
adjacent samples in the same spectral 
band and one sample each from the 
same location in each of three preced- 
ing spectral bands. Specifically, 


s(x-l,y,z)-s(x,y,z) 
s(x-l,y-l,z)-s(x,y,z) 
s(x,y - 1 ,z) - s(x,y,z ) 
k s(x,y,z- \)-s(x,y,z- 1) 
s(x,y,z - 2) - s(x,y,z - 2) 
s(x,y,z - 3) - s(x,y,z - 3) 


where s(x,y,z) is a mean value of pre- 
vious samples in the vicinity of x,y in 
spectral band z. The stream of e k val- 
ues calculated by use of this u* is fur- 
ther compressed by use of Golomb 
codes. 

In tests, the compression effective- 
ness of this algorithm was shown to 
be competitive with that of the best 
previously reported data-compres- 
sion algorithms of similar complex- 
ity. The table presents results from 
one series of tests performed on 
multispectral imagery acquired by 
NASA’s airborne visible/infrared 
imaging spectrometer (AVIRIS). 

This work was done by Matthew Klimesh 
of Caltech for NASA’s Jet Propulsion Lab- 
oratory. Further information is contained in 
a TSP (see page 1). 

The software used in this innovation is 
available for commercial licensing. Please 
contact Karina Edmonds of the California 
Institute of Technology at (626) 395-2322. 
Refer to NPO-4251 7. 
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Fast 


JPEG-LS 

Rice/USES 

Differential 

SLSQ 

SLSQ-OPT 

Scene 

Lossless 

ICER-3D 

(2-D) 

Multispectral 

JPEG-LS 

(Rizzo et al.)* 

(Rizzo et al.)* 


Cuprite 1 


4.89 

5.14 

7.13 

6.00 

5.44 

5.03 

4.90 

Cuprite 2 


5.02 

5.34 

7.50 

6.13 

5.58 

5.09 

4.97 

Cuprite 3 


4.92 

5.16 

7.16 

6.00 

5.45 

5.06 

4.92 

Cuprite 4 


4.98 

5.21 

7.16 

6.05 

5.51 

5.10 

4.96 

Jasper Ridge 1 


5.04 

5.41 

7.72 

6.17 

5.62 

5.06 

4.95 

Jasper Ridge 2 


5.02 

5.37 

7.67 

6.12 

5.59 

5.05 

4.94 

Jasper Ridge 3 


5.07 

5.47 

7.90 

6.19 

5.67 

5.10 

4.99 

Jasper Ridge 4 


5.07 

5.47 

7.87 

6.22 

5.67 

5.11 

5.00 

Jasper Ridge 5 


5.02 

5.39 

7.75 

6.14 

5.60 

5.06 

4.94 

Low Altitude 1 


5.37 

5.70 

7.81 

6.53 

5.97 

5.38 

5.30 

Low Altitude 2 


5.42 

5.76 

7.95 

6.58 

6.02 

5.40 

5.33 

Low Altitude 3 


5.30 

5.58 

7.57 

6.42 

5.88 

5.33 

5.23 

Low Altitude 4 


5.32 

5.58 

7.53 

6.42 

5.89 

5.37 

5.26 

Low Altitude 4 


5.37 

5.63 

7.60 

6.47 

5.91 

5.40 

5.30 

Low Altitude 6 


5.29 

5.56 

7.52 

6.42 

5.85 

5.34 

5.24 

Low Altitude 7 


5.29 

5.60 

7.64 

6.43 

5.88 

5.34 

5.24 

Lunar Lake 1 


4.99 

5.19 

6.98 

6.02 

5.49 

5.12 

4.99 

Lunar Lake 2 


4.94 

5.14 

6.96 

5.97 

5.44 

5.07 

4.93 

Moffett Field 1 


5.12 

5.48 

7.78 

6.24 

5.70 

5.15 

5.03 

Moffett Field 2 


5.11 

5.40 

7.57 

6.20 

5.60 

5.08 

4.98 

Moffett Field 3 


4.98 

5.12 

7.03 

5.96 

5.41 

4.96 

4.86 

Average 


5.12 

5.41 

7.51 

6.22 

5.68 

5.17 

5.06 


* F. Rizzo, B. Carpentieri, G. Motta, and J.A. Storer, "Low-complexity lossless compression of hyperspectral imagery via linear prediction," 
IEEE Signal Processing Letters, 1 2(2): 1 38-141, February 2005. 


Data From AVIRIS Images of various scenes were compressed by the present algorithm and by a number of other algorithms. The numerical entries are the 
numbers of bits per sample in the compressed data streams. 
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